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ABSTRACT 


We study the possibility of the existence of negative mass bubbles within a de Sitter space- 


time background with matter content corresponding to a perfect fluid. 


It is shown that there 


exist configurations of the perfect fluid, that everywhere satisfy the dominant energy condition, the 
Einstein equations and the equations of hydrostatic equilibrium, however asymptotically approach 
the exact solution of Schwarzschild-de Sitter space-time with a negative mass. 
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Introduction- The Schwarzschild solution [I] is an exact 
vacuum solution of the Einstein equations which contains 
one free parameter, the mass. The solution is singular at 
Schwarzschild coordinate r = 0. For a positive mass, the 
solution is called a black hole, and the singularity is hid- 
den behind a null surface called the event horizon. The 
singularity, rather than being concentrated at a point, 
actually occurs on a future space-like hypersurface. If 
an observer stays outside the event horizon, he is safe 
from the singularity. However, if he crosses over the null 
surface and enters the black hole he cannot avoid the sin- 
gularity, just as he cannot avoid aging, and in finite time 
he is ripped apart by the infinite gravitational stresses 
that occur at the singularity. On the other hand, for a 
negative mass, the singularity is naked as it has no event 
horizon cloaking it. However, it is point like, and oc- 
curs at a fixed spatial position. An observer can view 
the singularity from a distance, and in principle, every 
observer who chooses not to impinge on the singularity 
can avoid it forever. Indeed, the singularity of the nega- 
tive mass Schwarzschild solution appears to be relatively 
benign compared to the one of a positive mass. 


We are commonly exposed to singularities of the type 
exhibited by the negative mass Schwarzschild geometry. 
For example consider a point charge in electrodynamics. 
The corresponding electric field is in principle, singular 
at the location of the point charge, and the total energy 
is infinite. We are content with the understanding that 
from a large distance, a singular point charge is harmless, 
we can simply avoid it, and at close distance we expect 
that the singular nature of the charge will be smoothed 
out by a concentrated but non-singular charge density. 
The analogous situation in the context of general rela- 
tivity and the negative mass Schwarzschild solution, is 
to ask if there is a non-singular distribution of energy- 
momentum density that smooths out its singularity. This 
question was analyzed in [2] and it was concluded that for 
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asymptotically Minkowski space-time, it is actually easy 
to remove the singularity in the metric, but only at the 
expense of introducing other pathologies in the system, 
that are untenable. 


The physical definition of the mass is subjective. It can 
mean the total ADM mass [3] or the Bondi mass [4] or 
the Komar mass [5], depending on the symmetries and 
asymptotics of the metrics involved. We take the very 
physical view point of defining the mass as that parame- 
ter which would appear in the Newton dynamical equa- 
tions in the limit of slowly moving sources and weak grav- 
itational forces. This mass is exactly the Schwarzschild- 
like mass parameter that appears in the metric, as shown 
by Bondi in [6] and Luttinger [7]. 

A non-singular matter distribution that is physically 
reasonable and not pathological, must correspond to an 
energy momentum tensor that is non singular but addi- 
tionally, satisfies everywhere the dominant energy con- 
dition [8]. This condition is equivalent to the state- 
ment that no Lorentz observer can observe the energy- 
momentum to be moving out of the future directed light 
cone. It was observed in [2], that any non-singular mat- 
ter distribution, that asymptotically approaches the neg- 
ative mass Schwarzschild metric, must violate the dom- 
inant energy condition. The reason is because of the 
positive energy theorem [9], which states that in an 
asymptotically Minkowski space-time, any energy mo- 
mentum tensor that everywhere satisfies the dominant 
energy condition must give rise to a metric that has pos- 
itive ADM mass [3]. The negative mass Schwarzschild 
metric evidently has a negative ADM mass. Therefore, 
there cannot exist a smooth configuration of energy- 
momentum that everywhere satisfies the dominant en- 
ergy condition and asymptotically approaches the nega- 
tive mass Schwarzschild metric. Thus the negative mass 
Schwarzschild solution remains an unphysical solution 
and the question of its meaning is still unanswered. 


However, it is possible to evade the positive energy 
theorem by relaxing any of the assumptions that it re- 
quires. In [2], it was pointed out, that removing the con- 


dition that the space-time be asymptotically Minkowski 
is a useful reduction. In this letter we consider asymp- 
totically de Sitter space times. Positive energy theorems 
and the analysis of the positivity of the mass have been 
done in such space-times [IOHI4], but our analysis avoids 
their conclusions. The papers [10] prove the positivity of 
the mass associated with the conformal time translation 
Killing vector. Although positive, this mass is time de- 
pendent as far as the static time coordinate is concerned, 
and thus not what we would like to consider as the mass. 
The paper [II] and [I3] prove the positivity of the AD 
mass as defined by [12], but only if it is assumed to be 
positive on an initial surface. We do not make that as- 
sumption. The paper makes stringent conditions on 
the asymptotic space-time which also guarantees the pos- 
itivity of the mass, however we do not impose the same 
conditions. Finally the paper gives a nice counterex- 
ample to the conjecture that locally asymptotically flat 
space-times must also have positive mass. But none of 
these papers consider our specific, physical problem of a 
perfect fluid in an asymptotically de Sitter space-time. 

The Schwarzschild-de Sitter solutions are exact solu- 
tions of the Einstein equations in the presence of a con- 
stant, background energy-momentum density (or equiv- 
alently with cosmological constant), which correspond 
to an exponentially evolving space-time geometry with 
a singular, point like mass. Explicitly the metric is given 
by 
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The solutions contain two parameters, the value of the 
mass M, and the value of the constant energy density 
(or cosmological constant), A. Either can be positive or 
negative. Negative A gives rise to anti-de Sitter space- 
time, which, although of great interest recently [I6], al- 
ready on its own does not satisfy the dominant energy 
condition. Positive A gives rise to de Sitter space-time. 
The mass parameter, M, as expected describes a pos- 
itive or negative point-like mass depending on its sign. 
The positive mass is hidden behind an event horizon, 
while the negative mass singularity is naked, and occurs 
at a fixed spatial point. The positive energy theorem [9] 
is not applicable since the space-time is not asymptoti- 
cally Minkowski. The negative mass singularity can be 
smoothed out in a de Sitter background while maintain- 
ing the dominant energy condition everywhere, as was 
shown in [2], however, no attempt was made to find any 
type of energy-momentum which could give rise to the 
deformation that is required. In this paper we show that 
energy-momentum corresponding to a perfect fluid can 
be used to provide a suitable deformation. We find bub- 
ble like configurations which are non-singular, satisfy the 


dominant energy condition everywhere and give rise to a 
metric that asymptotically approaches the negative mass 
Schwarzschild-de Sitter geometry. 

Equations of motion- The metric of a spherically sym- 
metric space-time can be taken as 


ds? = B(r)dt? — A(r)dr? — r°dé? — r?sin?0d¢". (2) 


This metric is required to be a solution of the Einstein 
field equations: 
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where we take that c = G = 1. It has been commonly un- 
derstood that the left hand side of the equation Eqn. (3) 
depends on the geometry of the space-time and the right 
hand side depends on the matter content of the universe 
and its distribution. We wish to solve the Einstein field 
equations for matter corresponding to a perfect fluid, 
where the metric asymptotically approaches the exact so- 
lution corresponding to a negative mass Schwarzschild-de 
Sitter geometry. A perfect fluid is described by a stress 
energy tensor of the form 


Tuv = —P9pv + (p + p)U,U, (4) 


where p and p are respectively the pressure and the den- 
sity of the perfect fluid when U,, is its four vector ve- 
locity. We choose a frame of reference where the fluid 
is at rest, then this four vector velocity is defined by 
U, = (VB,0,0,0), so that it obeys U,U“ = 1. This 
means that the only non-zero terms of the stress-energy 
tensor are: 


T33 = sin? 0 T2 

(5) 
The Einstein equations in the presence of a perfect fluid 
are well known, for example in [I7], or in any other stan- 
dard reference on general relativity, however, we record 
them here for completeness. We get three independent 
equations: 
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These equations are underdetermined as they correspond 
to three equations among four fields in one independent 
variable, however, we will use them to determine the pres- 
sure p, the density p, the time component of the metric B 
as a function of the gravitational effective mass function 
m(r). A fourth equation corresponding to the condition 


of hydrostatic equilibrium is actually not an independent 
equation, and in fact a consequence of the covariant con- 
servation of the energy-momentum tensor, which in turn 
is equivalent to the consistency requirement of the Ein- 
stein equation with the Bianchi identities. It corresponds 
to: 
t 

a T 
p+p 
Usually, to render the system of equations determinis- 
tic, a further equation is introduced corresponding to the 
equation of state, p = p(p), which then picks out a unique 
solution. However we will not impose such an equation, 
we will instead parametrize 
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and choose m(r) explicitly, which has the interpretation 
as the mass inside the radius r. We then find the solution 
(numerically) for p(r), p(r) and B(r) dependent on our 
choice of m(r). It will be clear from our solution that 
an equation of state connecting p(r) and p(r) must exist 
implicitly. 

We can calculate the expression for the density through 


the operation Banh | za- | 2o H to obtain 
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and we can also get an equation for the time component 
of the metric B using the operation Han (6) ed td] 
to obtain 
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Finally, we can combine Eqn. (12) with Eqns. (6). (7). (8) 


to get 
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To be able to numerically solve the equations, we need 
to define the gravitational effective mass function m(r). 
We obtain a smooth function through the following steps. 
We start with the function 
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where x, y and a are parameters which are constrained 
by imposing that the interpolation be C? at x and y, and 
we have redefined M —> —M, see [18]. This implies 
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In the region r > y the space-time becomes the exact 
negative mass Schwarzschild-de Sitter space-time. Then 
we define a C% interpolation by 
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for a suitable value of w. 

Our choice for m(r) is motivated by the analysis done 
in [2] where it was shown that the dominant energy con- 
dition is satisfied if: 
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Although the analysis done in [2] does not exactly apply 
since there it was further assumed that B = 1—2m/(r)/r 
which we do not assume here, m(r) satisfies the condi- 
tions of Eqn. (17). The dominant energy condition for 
energy-momentum described by an ideal fluid simply re- 
quires 

Ar) 20 — p(r) = |p(r)| (18) 


which we will find are respected by our solution of the 
field equations. 
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FIG. 1: (color online) m(r) for the parameters z = .1, A = 
001, M = .01, y © 245. 


The de Sitter horizon is given by 
2m(r) —r =0 (19) 


which must be at a much larger radius than the size 
of the bubble given by y for a sensible configuration. 
The solution of Eqn. (19) must occur when r is in the 
Schwarzschild-de Sitter region of the space-time, r > y. 
Neglecting M then gives 


3 
AP yy—ra0are ye. (20) 


Comparing this with the expression for y from Eqn. (15) 


yields 
[6M [3 


Using these conditions combined with the expression 
of m(r), we get the numerical solution for the density 
and the pressure of the perfect fluid, p(r) and p(r), given 


respectively in Fig. 2} and Fig. (3). 
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FIG. 2: (color online) Density p(r) for the parameters x = .1 
G=1,A=.001, M = .01, y © 245. 
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FIG. 3: (color online) Pressure p(r) for the parameters x = .1 
G=1,A=.001, M = .01, y © 245. 


We note that the pressure is negative and interpolates 
from p(r) = 0 in the interior where the space-time is 
Minkowski to p(r) = —p(r), which is the normal solution 
for de Sitter space-time. The numerical solution for B(r) 
is not particularly simple to obtain because of numerical 
instability. Instead of directly numerically integrating 
Eqn. (8), or Eqn. (12), it is better to find an expression for 


AB , which comes straitforwardly from adding Eaa + 


Eqn. {7} . 

B: 
(In ABY = 87G (p(r) + p(r)) rA. (22) 
Numerically integrating Eqn. (22) yields the curve in 


Fig. (4). 

Thus our solution for the metric is given by Eqn.(10) 
for A(r), where m(r) is explicitly chosen as in Eqn. (16) 
and B(r) obtained from numerically integrating Eqn. (22) 
for (AB) and then trivially solving B(r) = e™4®)7A. 
To perform this integration, we need to first find p(r), 
which is given by Eqn. (11) and p(r), which is obtained by 
numerically integrating Eqn. (13) and shown in Fig. (3). 
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FIG. 4: (color online) In AB for the parameters x = .1 G = 1, 
A = .001, M = .01, y ~ 245. 


From Fig. (3) and Fig.(2) we can see the combination 
(p(r) + p(r)) vanishes as r becomes large. Thus the ex- 
pectation is that AB becomes a constant for large r as 
can be seen from Fig.(4). This constant can be taken 
to be 1, although it is not so in Fig.) as the integra- 
tion of Eqn. (22) is only determined up to a constant. 
Notice that the Einstein equations, Eqns. (6) ,(7), (8) .(9) 
all are homogeneous in B. If B is a solution then so is 
aB for any constant (positive) a. Thus asymptotically, 
B = 1/A and the solution approaches the exact negative 
mass Schwarzschild-de Sitter solution. 

If we plot p(r)/p(r) we find the curve given by Fig. (5), 
which is bounded between 0 and -1. It is clear then, 
within the limits of numerical accuracy, that the con- 
straints of the dominant energy condition, Eqn. (18) are 
satisfied. 
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FIG. 5: (color online) p(r)/p(r) for the parameters x = .1 
G =1, A = .001, M = .01, y © 245. 


Conclusions- We have shown that there exist very 
physical configurations of an ideal fluid which give rise 
to solutions of the Einstein equations that correspond 
asymptotically to negative mass Schwarzschild-de Sitter 
space times. The energy-momentum tensor that gives 
rise to such space times is perfectly physical, it every- 
where satisfies the dominant energy condition. Since the 
space time is not asymptotically flat, we evade the posi- 


tive energy theorems which would not allow for negative 
mass. Negative mass configurations therefore can exist in 
de Sitter backgrounds, exactly as have been proposed for 
the inflationary phase of the early universe. If a mecha- 
nism for production of pairs of particles with positive and 
negative mass can be determined, in the early universe 
there would be a plasma of positive and negative mass 
particles. Such a plasma would in principle cause an ef- 
fective screening of gravitational waves, being essentially 
opaque for frequencies below the plasma frequency. 
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